Dynamics, dynamic soft elasticity and rheology of smectic-C elastomers 
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We present a theory for the low-frequency, long-wavelength dynamics of soft smectic-C elastomers 
with locked-in smectic layers. Our theory, which goes beyond pure hydrodynamics, predicts a 
dynamic soft elasticity of these elastomers and allows us to calculate the storage and loss moduli 
relevant for rheology experiments as well as the mode structure. 
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Smectic elastomers [l| are rubbery materials that have 
the macroscopic symmetry properties of smectic liquid 
crystals 2.J. They are sure to have intriguing properties, 
some of which have already been studied experimentally 
and/or theoretically Very recently, seminal progress 
has been made on smectic-C (SmC) elastomers forming 
spontaneously from a smectic- A (SmA) phase upon cool- 
ing. Hiraoka et al. Q for the first time produced a mon- 
odomain sample of such a material and carried out ex- 
periments demonstrating its spontaneous and reversible 
deformation in a heating an cooling process. Also very 
recently, it was discovered theoretically that such a ma- 
terial exhibits the fascinating phenomenon of soft elastic- 
ity 0, i.e., certain elastic moduli vanish as a consequence 
of the spontaneous symmetry breaking wherefore strains 
along specific symmetry direction cost no elastic energy 
and thus cause no restoring forces. 

On one hand, due to the aforementioned experimen- 
tal advances, dynamical experiments on soft SmC elas- 
tomers, such as rheology experiments of storage and loss 
moduli or Brillouin scattering measurements of sound ve- 
locities, seem within reach. On the other hand, there 
exists, to our knowledge, no dynamical theory that 
could be helpful in interpreting these kinds of experi- 
ments. Here we present a theory for the low- frequency, 
long- wavelength dynamics of soft SmC elastomers with 
locked-in layers that goes beyond pure hydrodynamics. 
As in standard elastic media and nematic elastomers 
a purely hydrodynamical theory of SmC elastomers in- 
volves only a displacement field u and not the Frank di- 
rector n, which relaxes to the local strain in a nonhydro- 
dynamic time r„. We go beyond hydrodynamics, by in- 
cluding n in our theory, because dynamical experiments, 
like rheology measurements, typically probe a wide range 
of frequencies that extends from hydrodynamic regime to 
frequencies well above it. 

Smectic elastomers are, like any elastomers, perma- 
nently crosslinked amorphous solids whose static elastic- 
ity is most easily described in Lagrangian coordinates in 
which x labels a mass point in the undeformed (refer- 
ence) material and R(x) = x + u(x), where u(x) is the 



displacement variable, labels the position of the mass 
point x in the deformed (target) material. Lagrangian 
elastic energies are formulated in terms of the strain ten- 
sor u which, in its linearized form, has the components 

u ij = hiVij + Vji)> where rjij = djUi are the components 
of the displacement gradient tensor rj. 

The elastic energy density / of the SmC elastomers of 
interest here can be divided into two parts y| 

/ = /u + /u,„, (1) 

where / u depends only on u and / u , n describes the depen- 
dence of / on the Frank director n including its coupling 
to the displacement variable, fn the following we choose 
the coordinate system so that the z-axis is parallel to the 
director of the initial SmA phase and the x-axis is par- 
allel to the direction of tilt in the resulting SmC phase 
so that the equilibrium director characterizing the unde- 
formed SmC phase is of the form n° = (c, 0, yl — <?) 
with c being the order parameter of the transition. With 
these conventions, / u can be written in the same form 
as the elastic energy density of conventional monoclinic 
solids 0, 

/u = 2 Cxyxy U xy + C xyzy U xy U zy + ^ C zyzy uj, y 
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but with constraints relating the three elastic constants 
in the first row. These latter constants can be expressed 
in terms of an overall elastic constant C and an angle 
9, which depends on the order parameter c, as C xyxy = 
C cos 2 8, Cxyzy — C cos 9 sin 9, and C zyzy = C sin 2 9. Ne- 
glecting contributions from the Frank energy, / U) „ can be 
stated as 

/ u ,n = \ A [Q y + au xy + l3uy X } 2 (3) 

where A is a coupling constant and where a and (3 are 
dimensionless parameters. The variable Q y stands for 
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Qy 



cr/Ayx - vl 
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where rjAyx and rj Ayz 



are components of the antisymmetric part of r\. 

Equations (JTJ to © imply that SmC elastomers are 
soft under static conditions. Imagine that Q y has re- 



laxed locally to Q y 



f3u yx so that / is ef- 



fectively reduced to / u . Then, due to the above rela- 
tions among the elastic constants, deformations charac- 
terized by Fourier transformed displacements u || e y and 
wavevectors q || ea = (— sin 9, 0, cos 9) or, alternatively, 
u || ea and q || e y cost no elastic energy and hence cause 
no restoring forces. The effects of soft elasticity are more 
evident in a coordinate systems rotated though 9 about 

i y , and 
vanish 



the y-axis in which e x > — (cos 9, 0, sin 9), 



■y 



e z i = f»2. In this system, C x >y'z'y' and Cz' y i z 'y' 
impyling there is no energy cost for shears u z ' y > in the 
z'y'-plane. If the elastomer is crosslinked in the SmC 
phase, these moduli become nonzero. Here, we will not 
consider such semi-soft SmC elastomers. 

Now, let us formulate our dynamical theory. Dynami- 
cal equations for u and n can be derived using standard 
Poisson-bracket approaches Q, with the result 

STi 

hi = Xijk djiik-T — , (4a) 

, „ STi STi _ „ . , . 

pui = \ kji dj- h Uij M djdiUk , (4b) 

0Tl k OUi 

where Ti is the elastic energy of the system, v^ju is the 
viscosity tensor, T is a dissipative coefficient with dimen- 
sions of an inverse viscosity, and 

X ijk = \ A (Sj- n k + Sf k n k ) + \ (Sfi n k - 8j k n k ) , (5) 

where Sfj — §ij — riiiij. As they stand, Eqs. J3J are valid 
for any liquid crystal elastomer with a defined director, 
like e.g., nematic, SmA and SmC elastomers. To describe 
SmC elastomers we have have to specify Ti and Vij k i ac- 
cordingly. From the above it is clear that Ti = J cPxf 
with / as given in Eq. IjTJl. The viscosity tensor entering 
here is that of a monoclinic system. It has 13 indepen- 
dent components and it can be parameterized, as we do, 
so that the entropy production density Ts takes on the 
same form as Eq. (J3J) with the elastic constants C^u re- 
placed by viscosities Vijki and with tiy replaced by iiij. 
The Viju depend on the order parameter c. v xxxz , v yyxz , 
and v zzxz vanish at the SmC to Sm^4 transition and are 
therefore expected being smaller than the remaining vis- 
cosities for c small Q- 

A smectic elastomer is characterized in general by re- 
laxation times associated with director relaxation and 
with other modes, which we will simply refer to as elas- 
tomer modes. For frequencies w such that lute -C 1, 
where te is the longest elastomer time, the viscosities and 
r are practically frequency independent. When lute > 1, 
however, the viscosities v^u and V develop a non-trivial 
frequency dependence. In the following we will consider 
in detail only the case r„ 3> te and lute *C 1. 



As mentioned above, Eq. Q omits contributions from 
the Frank elastic energy for director distortions, which 
are higher order in derivatives than those arising from 
network elasticity. Without the Frank energy, our dy- 
namical theory misses diffusive modes along certain sym- 
metry directions where sound velocities vanish. Including 
the Frank energy, on the other hand, makes the equations 
of motion considerably more complicated. To keep our 
presentation as simple as possible, we will therefore, for 
the most part, exclude the Frank energy. When it comes 
to stating results for the aforementioned diffusive modes, 
however, we will include Frank contributions in order to 
present complete results. 

From Eq. H4a[) we can derive an equation of motion for 
Q y . In frequency space, this equation can be written as 



Qy — 



1 + ILUT 3 
1 — ILUT\ 



1 + iLUT 2 
1 — ILUT\ 



(0) 



where we have introduced the relaxation times t\ = 
1/(TA), t 2 = XVT^/iTAp), t 3 = Xc/(TAa). As we 
will see further below, our dynamical equations predict 
nonhydrodynamic modes with a decay time ( "mass" ) t\ 
which implies t\ = r n . 

With help of Eqs. ® and (gEJl we derive effective equa- 
tions of motions in terms of the displacements only which 
can be cast as 

pLu 2 Ui = -djdijiuj) (7) 
with a symmetric stress tensor a given by 



a^uj) = C^(lu) u u + C^ ixz (uj) u xz , 
(J xz (lu) = \ C xzv (lu) u y + \ C xzxz (lu) u xz , 

<x«M = |cf» 



(8a) 
(8b) 
(8c) 



where we use a compact notation with indices /i, v run- 
ning over xx, yy and zz and indices £, x running over 
xy and zy. Cij k i (lu) with no superscript R stands for 
Ciju(Lu) = C'ij k i — iujVijki- The superscript R indicates 
that certain elastic moduli are renormalized by the re- 
laxation of the director. These are 



Cf x {w) = C £x -iuJV £x 



C 



ILUT\ 

1 — ILUT\ 
,2\ 



■ AA 



ix 



OK), 



(9) 



r _1 A ?x , and 



with renormalized viscosities = v £x 
where A xyxy = a 2 (l + r 3 /Ti) 2 , A xyzy = af3(l+T 3 /n)(l + 
T 2 /n) and A zyzy = /3 2 (1 + t 2 /ti) 2 . 

The frequency dependence of the elastic moduli in 
Eq. JSJ can be determined, in principle, by rheology mea- 
surements of the corresponding storage and loss moduli. 
The unrenormalizcd moduli Cij k i(uj) lead to conventional 
storage and loss moduli C' ljkl — C t j k i and C" jkl = Luu ljk i 
that are, as in conventional rubbers, respectively con- 
stant and proportional to lu at low frequencies. The 
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FIG. 1: Log- log plot of the reduced storage and loss moduli 
C' ix [oj)/C tx (solid line) and C^(w)/Cf x (dashed line) versus 
the reduced frequency u>T n as given respectively by the real 
and negative imaginary parts of Eq. © for r n S> te- For the 
purpose of illustration we have set, by and large arbitrarily, 
v (x /(T n C ix ) = 1 and AA ix /C\ x = 10 2 . 

rcnormalizcd moduli, on the other hand, have the poten- 
tial for much more interesting rheology behavior. One 
consequence of Eq. © is that SmC elastomers could ex- 
hibit so-called dynamic soft elasticity 9]. To highlight 
this phenomenon, let us switch briefly to the rotated ref- 
erence space coordinates x',y',z'. Then, Eq. 10 implies 
that 

Cy'z'y'z' M = ~ iuJ Vy'z>y>z> ~ ^ AA y > z ,y, z , (10) 

with Vy' Z 'y' z ' = sin 2 6v xyxy - s\n29v xyyz + cos 2 9v yzyz 
and an analogous expression for A y i z i y t z i . C^; , , z , (oj) 
is of the same form as Eq. {Hfl) . Thus, Cy, z , y , z ,{oS) and 

^x'y'y'z' ( w ) van ish m the limit u> — > where we recover 
true soft elasticity. At non-vanishing frequency the sys- 
tem cannot be ideally soft but it can be nearly so for lu 
small. This type of behavior was first predicted for ne- 
matic elastomers, where it has been termed dynamic soft 
elasticity || . The storage moduli for u xy and u yz strains, 
Eq. © , in the original coordinate system are nonzero at 
zero frequency. Their behavior for ui > 0, like that of 
semi-soft nematic elastomers 0, 0, depends on t u /te- 
If r n 3> te, the storage moduli exhibit a step and the 
corresponding loss moduli an associated peak at wr„ ~ l 
as shown in Fig. 2] if r„ ~ te, or t„ < te, this is not the 
case. The storage moduli C' y , z , y , z ,(uj) and C' x , y/y , z ,(uj) in 
the rotated frame are zero at zero frequency. In a semi- 
soft SmC, however, they will exhibit behavior similar to 
that Fig. H for r„ 3> te- I n nematic elastomers, there is 
still some controversy [lfj about whether the r„ te 
regime has actually been observed in experiments. It 
would be interesting to see if it might exist in SmC elas- 
tomers, where te might be shorter than it is in nematics 
because of the smectic layers. 

To assess the mode structure of SmC elastomers, we 
start with an analysis of propagating sound modes in the 
dissipationless limit. The sound modes have frequencies 
w(q) = C("d,(p)q, where q = |q| and where z? and ip are 
the azimuthal and polar angles of q in spherical coordi- 




FIG. 2: Schematic plots of sound velocities: (a) spherical plot 
of mode pair (i) only; (b) and (c) polar plots in the y'z'- and 
x'z'- planes, respectively, of all 3 sound-mode pairs. 

nates. Their sound velocities C($, <p), as calculated from 
Eq. 10 with the viscosities Vijki set to zero, are depicted 
in Fig. El There are 3 pairs of sound modes. One of 
these pairs (i) is associated with the soft deformations 
discussed above. Its velocity vanishes for q along e y and 
i-2 so that when viewed in the y'z'-plane it has a clover 
leave like shape. The remaining 2 pairs are associated 
with non-soft deformations. In the incompressible limit, 
these pairs become purely transverse (ii) and longitudinal 
(iii), respectively. In the y'z' and x'z'-planes, their veloc- 
ities are non-vanishing in all directions. Note that, since 
the velocity of pair (i) vanishes in directions where the ve- 
locities of the other modes remain finite, SmC elastomers 
are, like nematic elasomers 0, potential candidates for 
applications in acoustic polarizers. 

Having found the general sound-mode structure in the 
nondissipative limit, we now turn to the full mode struc- 
ture in the incompressible limit. In the softness-related 
symmetry directions the modes are of the following types: 

(i) non-hydrodynamic modes with frequencies 

W m = -iTi 1 + iD m q 2 (11) 

with a zero-q decay time n and diffusion constants D m , 

(ii) propagating modes with frequencies 

u p = ±Cq~iD p q 2 (12) 

with sound velocities C and diffusion constants D p , and 

(iii) diffusive modes with frequencies 

uj d = -W d q 2 ± ^-(D d q^ + Bq* (13) 
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with diffusion constants D d and bending terms B that are 
missed if the Frank energy is neg lected. For B/D 2 d <C 1 
the diffusive modes split up into slow and fast modes 



where i>3„ 



is the renormalized version of v, 



and 



uj s = -iB/{2D d )q 2 



l2D d q 2 



(14) 



Specifics of the sound velocities, the diffusion constants 
and the bending terms are given in the following. 

First, let us consider the case that q lies in the xz- 
plane. In this case the equation of motion for u y decou- 
ples from the equations of motion for u x and u z . This 
equation produces a set of transverse modes with u || e y . 
There is a non-hydrodynamic mode with 



D m ,y = (4p) 1 \yfZ-. 



I <fa + \J v zyzy v zyzy Qz 



(15) 



where % = qi/q, and there are propagating modes with 

C y = ^C/{A P )\cos6q x +sm6q z \ = C / {Ap)\q x ,\ , 

(16a) 

Dp.y = (Sp)' 1 [y% xv ol + 2^ yzy q x q z + v? vxy fi] ■ (16b) 



In the soft direction, i. e. for q 



these propagating 



'y'zy ~ — ~ " — ^z'y'z'y' 

where the K's are once more bending moduli depending 
on the Frank constants, c, A, a, and j3. Finally, there is 
a pair of propagating modes polarized along e x * with 



+ [C xxxz - C zzxz ]q x q z (q 2 z - q 2 x ) + \C xzxz {q 2 z - q 2 x )} 
D p .t = (2p) _1 [[v xxxx + v zzzz - 2v xxzz ]q 2 x q 2 z (18b) 
+ \Vxxxz - v zzxz ]q x q z (q 2 z - q 2 x ) + \v xzxz {q 2 z - q 2 x )] , 

and u || e T where e T = (q z ,0, -q x ). 

Finally, we turn to the case q || 
of longitudinal propagating modes with u 
suppressed in the incompressible limit. There is a non- 
hydrodynamic mode with 



There is a pair 
, that is 



D m , xz = (ip)- 1 [vl 



+ VV„*„ - Vz 



(19) 



where u lays in the xz-plane with u x = a{r\ + ts)/[P(t\ + 
T2)] u z . There is a pair of elastically soft diffusive modes 
with polarization u || e z i with 

D d , z , = (Sp)- 1 v%, zlyl , (20a) 
B z > = p' 1 [sin 2 8K 6 - sin 29 K 8 + cos 2 9 K 7 ] , (20b) 



C xl =^C/{Ap) : 



D p , x , = (Sp)- 1 uX'x'y' . 



(21a) 
(21b) 



with V x>y'x<y 



cos" 



sin 29 v* yzy 



sin 2 9 v?„. 



xyxy • ""^ " xyzy 1 zyzy 

In summary, we have presented a theory for the low- 
frequency, long-wavelength dynamics of soft SmC elas- 
tomers. This theory predicts that, at least in an idealized 
limit, SmC elastomers possess the fascinating property of 
dynamic soft elasticity. Though the equations of motion 
are complicated, the resulting mode structure is with re- 
spect to the softness related symmetry directions nicely 
symmetric and it has, when visualized, a certain beauty. 
We have calculated various dynamical quantities, such as 
storage and loss moduli and sound velocities that are, in 
principle, accessible by experiments and we hope, that 
our theory encourages such experimental work. 

Support by the National Science Foundation under 
grant DMR 0404670 (TCL) is gratefully acknowledged. 



modes become diffusive with D d y — D p y and 

B y = p- 1 [K ig 4 x + K 2 q 4 z + K 3 f x f z + 2K iq 3 x q z + 2K 5 q x q 3 z ] , 

(17) 

where the J^'s are bending moduli that are combinations 
of the usual Frank elastic constants, the order parameter 
c as well as A, a, and (3. The equations of motion for 
u x and u z can be solved by decomposing (u x ,u z ) into a 
longitudinal part m along q and a transversal part ut- 
In the incompressible limit ui vanishes. The equation of 
motion for u T produces a pair of propagating modes with 

C T = VV~ P { [C xxxx + C zzzz - 2C xxzz ]q 2 x f z (18a) 
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